Continuum in the spin excitation spectrum of a Haldane chain, observed by neutron 

scattering in CsNiCl3. 
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The spin excitation continuum, expected to dominate the low-energy fluctuation spectrum in the 
Haldane spin chain around the Brillouin zone center, q = 0, is directly observed by inelastic magnetic 
neutron scattering in the 5* = 1 quasi-lD antiferromagnet CsNiCl3. We find that the single mode 
approximation fails, and that a finite energy width appears in the dynamic correlation function 
S(q,u>) for q < 0.57T. The width increases with decreasing q, while S(q,u>) acquires an asymmetric 
shape qualitatively similar to that predicted for the 2-magnon continuum in the nonlinear cr-model. 
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PACS numbers: 75.10.Jm, 75.40. Gb, 75.50.Ee 

Experiments on low-dimensional magnets are an in- 
valuable source of insight into the fundamental proper- 
ties of quantum spin models and field theories, many of 
which are still awaiting theoretical solutions. A vivid ex- 
ample is provided by the one-dimensional (ID) 5 = 1 
Heisenberg antiferromagnet (HAFM). It was conjectured 
by Haldane [jlj , and then confirmed by experiments |^-|^] 
and numerical studies P~pl]], that its ground state (GS) 
exhibits a robust quantum disorder, characteristic of the 
ID nonlinear a model (NLctM). In striking contrast with 
the almost ordered "marginal liquid" state of the S = 1/2 
chain |l^] , the Haldane GS is a " dense spin gas" , with a 
short correlation length of £ w 6 lattice repeats, and a 
gap Ah w 0.41J (J is the exchange coupling) in the spin 
excitation spectrum. The excitations' spectral weight is 
concentrated in a long-lived massive triplet mode in the 
vicinity of the Brillouin zone (BZ) boundary q = 7r, some- 
what similar to the spin-wave doublet in the classical 
1 AFM with easy-axis anisotropy Any remains of 
the spectacular continuum at q w it, found in the S — 1/2 
HAFM chain is attributed to 3(5,7,...)-magnon 

states, and predicted to be extremely weak ^,0,0. Al- 
though no quantitative theory of the excitation spectrum 
in the whole BZ exists, a continuum is expected in the 
vicinity of the zone center q — 0, where the lowest energy 
excitations are pairs of q ~ n magnons. 

The two-magnon nature of the lowest excited state in 
the ID S = 1 HAFM at q w was first stipulated by 
the Monte Carlo (MC) numerical study ||. It estab- 
lished, that the gap in the magnon spectrum at q = is 
about twice as large as that at q = ir. It was later sup- 
ported by MC simulations on larger systems and 
exact diagonalization (ED) for rings with up to 18 spins 
||[)]]. A semi-quantitative illustration of the lattice pe- 
riodic quasiparticle dispersion, crossing the two-particle 
energy at a q < 7r/2, was obtained in the physically mo- 
tivated variational theory [051, based on the mapping of 



the spin Hamiltonian on a system of spinless fermions. 
Not only does this simple theory correctly capture the 
fermion-like nature of the q m -k magnons observed in ex- 
periment pi , but it also gives a correct functional form 
for the single-magnon dispersion, 



e(q) = JA^+w 2 sinV 
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which for q > 7r/2 is in a very good agreement with 
MC ||. In fact, the dispersion, measured by neutron 
scattering in the Haldane chain compound NENP ||, 
was found to coincide with ([!]) in a large part of the 
BZ, at q > 0.37T, if the parameter a is adjusted from 
a w 2.5 J Jl5j to a w 1.45 J. In this paper we use ([!]) 
to describe the ID excitation spectrum in CsNiChj (solid 
lines in Fig. 1,3). We find that it gives an excellent fit to 
the measured excitation energies, bringing an end to the 
existing controversy [Q. Using J = 2.275 meV, indepen- 
dently determined from the high field magnetization p7[ , 
we obtain v — 2.49(4)J, in a remarkable agreement with 
the numerical result |11| . Most importantly, we find that 
spectral density deviates from the SMA at q < 0.57T, and 
present the first resolved measurement of the continuum 
part of the excitation spectrum of a Haldane spin chain. 

In spite of extensive theoretical evidence for exis- 
tence of the q sa continuum in the spectrum of a 
ID S = 1 HAFM, its limits and extent remain un- 
clear and controversial, while its experimental charac- 
terization is considered practically impossible because 
of the rapid decrease of the static spin structure factor 
S(q) = J S(q,oj)d(fiLij) at small q. This follows from the 
first moment sum rule [i~8| , which establishes an exact 
relation between S(q) and the excitation average energy 
HQ)) = J(hLu)S{q,Lu)d(huj)/ J S(q,uj)d(hLu). For the ID 
HAFM, 
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where E GS = -1.40(0) J |l0| is the GS energy per site. 
For a gapful spectrum S(q) of (0) vanishes ~ q 2 at q — > 0. 
Where the SMA holds, 5(g) is uniquely determined by 
the dispersion, (e(q)) = e(q). In addition, experimental 
characterization of the continuum is challenged if the ex- 
citation spectrum is split by the anisotropy, as in NENP 
and related Ni-organic model compounds, or if the exci- 
tations acquire a temperature-activated damping. 

GsNiCla is one of the most isotropic and best studied 
quasi-lD S = 1 HAFM model compounds (see |^,|| and 
references therein). It has a hexagonal crystal structure, 
space group P63/mmc; at T = 1.5 K the lattice spac- 
mgs are a = 7.12A, c = 5.9A. Chains of Ni 2 + 10ns run 
along the c axis and form a triangular lattice in the a — b 
plane. There are two equivalent ions per c spacing, so 
that Q — (h, k, I) in reciprocal lattice units (rlu) corre- 
sponds to q — ttI in the ID BZ of a chain. A very reliable 
estimate for the in-chain exchange coupling, J = 2.275 
meV, is obtained from the measured spin-flip (saturation) 
field H s = AJgfi B S = 73.5 T @. A spin-flop (reorien- 
tation) field H s f ss 1.9 T implies a negligible single- ion 
anisotropy, D w 0.002 J. 

We measured the dependence of the energy spectrum 
of the spin dynamic structure factor on the wavevector 
transfer along the chains for a CsNiCFj sample made of 
two large crystals with total mass 6.4(1) g, at T = 1.5 
K. The sample was mounted with the (h, h, I) zone in the 
scattering plane, and had effective mosaic spread < 1°. 
To increase the data collection rate, the SPINS 3-axis 
cold neutron spectrometer at NIST Center for Neutron 
Research was equipped with a position sensitive detec- 
tor (PSD), matched in size to the large PG(002) flat 
analysing crystal. The analyser's central energy was fixed 
at E^ = 4.2 meV; its angular acceptance was w 9°. Cal- 
ibration of the neutron final energy and the sensitivity 
across the PSD was done using elastic incoherent scat- 
tering from Vanadium. Beam divergence was defined by 
the 58 Ni neutron guide and a 80' radial collimator in front 
of the PSD. To perform energy scans at constant q = nl, 
the wavevector transfer along [110] was varied, so that 
chains were always aligned with the analyser [002] Bragg 
wavevector ta- Typical variation of h, imposed by this 
condition, is illustrated by the right and top axes of Fig. 
1, which shows the contour plot of the raw spectral den- 
sity of the scattering intensity (~ I(q,E)/ J I(q, E)dE), 
with the g-dependent flat background subtracted. It is 
evident from Fig. 1, that the spectrum acquires a finite 
width in energy at / < 0.5. Although the effect is some- 
what exaggerated by the resolution, careful accounting 
for the latter shows that intrinsic width, where it is non- 
zero, accounts for ~ 2/3 of the total width. Curves show 
the prediction of Jl5| , obtained from (|j) with J = 2.275 
meV, A H = 0.41J , v = 2.49 J, and a = v. 

For the quantitative analysis of the energy dependence 
of the measured cross-section, we use the normalized 



scattering function of the damped harmonic oscillator 
(DHO), parametrized in terms of the position ujq, and 
the full width at half maximum (FWHM) of the corre- 
sponding antisymmetrized Lorentzian peak Alterna- 
tively, to describe the asymmetric peak shape at lower 
q, we use a " half-Lorentzian" truncated DHO (TDHO), 
which is obtained by multiplying the above by 9(lu 2 — ujq). 
The TDHO fits to several constant-g scans are shown in 
Fig. 2. As quantified by X 2 in Fig. 3(d), TDHO is in 
better agreement with experiment at q < 0.57T, where a 
sharper onset of the scattering at low energies, character- 
istic of a continuum, is observed. Note, that an opposite 
asymmetry Q would result from the interplay of the dis- 
persion with the instrument resolution. 
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FIG. 1. Contour plot of the measured spectral density of 
magnetic scattering, reconstructed from 9 constant-q scans 
via linear interpolation. Scale on the right shows variation 
with energy of the wavevector transfer perpendicular to the 
chain at / = 0.5, scale on the top - its variation with / at 
E — 3 meV. Ellipses are the half maximum contours of the 
instrument resolution function, calculated at E — 3 meV. 
Solid curve is the single-magnon dispersion (0), dashed line 
shows the lowest energy of the two non-interacting magnons 
with given total q — irl, dotted line is e(q) = using. 

Peak parameters, refined for both DHO and TDHO 
fits, are detailed in Fig. 3. The dispersion of the cen- 
ter of mass of the excitation spectrum, captured by the 
position of the DHO peak (circles in Fig 1(a)), is nicely 
described by ([!]). For J = 2.275 meV, the best fit, shown 
in the figure, gives A H = 0.34(6)J, v = 2.49(4) J, and 
a = 1.1(4). While at q < 0.9tt the spectrum is not very 
sensitive to Ah, the agreement of the magnon velocity 
v with calculations Jll| is impressive. The energy in- 
tegrated intensity of both DHO and TDHO peaks, Fig. 
3(b), is in good agreement with the sum rule (En. A slight 
deviation, concomitant with non-zero peak energy width 
(Fig. 3(c)), is due to the extra intensity in the unphys- 
ical high-energy tail, and can be completely removed by 
imposing an upper boundary on the continuum. 
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FIG. 2. Energy spectrum of the inelastic scattering, mea- 
sured in CsNiCl3, for several q — Qc/2 — irl. The wavevector 
component in the a — b plane q± = Qa + Qb — irhy/3 is de- 
fined by the scattering geometry c || ta (Fig. 1). Solid circles 
show the non-magnetic background. Thick lines are the reso- 
lution-corrected fits to the cross-section, based on the TDHO 
scattering function. Shaded peaks are the same cross-section, 
multiplied by the resolution volume, and illustrate the "de- 
convoluted" intensity. Dotted lines show the cross-section for 
the 2-magnon continuum in the NL<r- model , corrected for 
the non-relativistic dispersion ([j]). 

The peaks at q < 0.57T, are very similar in shape and 
intensity to the two-magnon continuum, calculated for 
the ID NLcrM @ , which correctly describes the S = 1 
HAFM chain in the vicinity of q = it. We show the NLctM 
prediction by dotted lines in Fig 2. Even though we used 
the realistic magnon d ispersion (p|) , instead of the rela- 
tivistic form e(q) = y/ + (vq) 2 , q — qmodTr, which 
fails for q > 0.27T, the calculated continuum is still too 
high in energy. In fact, this discrepancy does not simply 
show the limitation of the cr-model, but also reveals the 
importance of magnon interactions. Indeed, a careful ex- 
amination of Figs. 1, 3(a) shows, that the lower boundary 
of the observed continuum lies below the lowest energy of 
two non-interacting magnons. This could be understood, 
if there is an attraction between the quasiparticles, which 
form the continuum states. 

To understand how our findings compare with existing 
results, it is important to realize that the extent of the 



observed continuum is rather small, as it should be, ac- 
cording to the calculations Q. At q = O.Stt, the FWHM 
is w 0.3 J, or only w 15% of the peak energy, and re- 
quires high resolution and low temperature, T <C 0.3J, 
to be detected. Previous attempts to measure the in- 
chain dispersion in CSN1CI3 [0 were, on the contrary, 
done with rather coarse resolution of the thermal neu- 
tron spectrometer, and at T ~ 0.4J, where temperature 
damping dominates the spectral width. It is more inter- 
esting to compare our data with the neutron scattering 
study in NENP |3| , where the authors measured the dis- 
persion for q > 0.37T at T < 0.01J, and found no appre- 
ciable deviation from the SMA. Consequently, a picture 
for the excitation spectrum of the ID S = 1 HAFM, 
where a single-mode dispersion merges into a broad, but 
unmeasurable, continuum at q < 0.37T, became broadly 
popular ^,0. The results of || are, however, easily 
reconciled with our data, if we note that, to observe 
magnetic scattering at small q, the authors had to re- 
lax the spectrometer resolution so much that almost a 
quarter of the ID dispersion was within their wavevector 
acceptance. The effect we observe is simply unresolved 
in such a measurement. The splitting of the spectrum by 
single-ion anisotropy, which is rather large in NENP [jl6| , 
presented another obstacle to the characterization of the 
continuum. A crossover to a continuum at q « 7r/2 is 
also supported by numerical calculations |7|, |ll| , |l(| . 

Finally, the main disadvantage of CsNiCl3 as a model 
ID HAFM is the supercritical inter-chain coupling J±, 
which leads to a 3D order with the propagation vector 
Q Q = (1/3, 1/3, 1), at T N fa 4.8 K. The relevance of the 
Haldane conjecture in this situation was the subject of 
a long-lasting controversy. It was established that a gap 
opens in the spin excitation spectrum at T > TV, and 
in this so-called " ID phase" it recovers the features of a 
Haldane chain. However, since TV ~ 0.2J, the spectra 
also acquire significant temperature broadening H. In 
fact, we argue that 3D order in CsNiCla is so weak, that 
even at T « it causes no significant change in the spec- 
trum throughout the better part of the BZ. This agrees 
with the simple physical argument that excitations whose 
energies are sufficiently large compared to J± are not sen- 
sitive to it. Indeed, elastic Bragg intensity, corresponding 
to the ordered spin value (S) = 0.5, observed in experi- 
ment |(|, accounts for only (S) 2 /{S{S + 1)) = 12.5% of 
the spin fluctuation spectrum. In the ID S — 1 HAFM 
this fraction of the spectral weight is concentrated in the 
tiny region of the BZ, at \q\ > 0.987T. Experiments also 
show, that the Haldane gap triplet mode is split in the 
3D ordered phase only in the close vicinity of the mag- 
netic Bragg peaks, at \q\ > 0.9-7T, where the Goldstonc 
acoustic magnons appear 

To quantify our arguments, we consider corrections to 
the S{q) of (g) for the ordered quasi-lD HAFM. The 
right-hand side is modified in three ways. First, the ex- 
change energy of an individual chain changes on account 
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FIG. 3. Wavevector dependence of the excitation spectrum 
in CsNiCl3. Full and open circles show the parameters ob- 
tained for the asymmetric TDHO and symmetic DHO line- 
shapes, respectively. Crosses are ED results Q. (a) Shaded 
bars start at the lower boundary of the asymmetric TDHO 
peak, and extend by its FWHM. Circles show the Lorentzian 
peak position of the symmetric DHO fits. Lines are the same 
as in Fig. 1, but with Ah = 0.34J, a = 1.2,1. (b) Static 
structure factor. Line is the SMA result for the dispersion 
shown in (a), (c) Peak (continuum) energy width, (d) Stan- 
dard square deviations for the TDHO and DHO fits. 

of the static order. In the mean field (MF) approxima- 
tion, at (S) < 0.5, it increases by (S) 2 /(2xn), which is 
less than 0.5%, due to the large static staggered suscep- 
tibility of a Haldane chain Xn ~ 22/ J [ppl|]. Second, the 
numerator has to be amended by adding the interchain 

correlation, J2p, T± J ^( X ~ cos(q ± • T±))(S% L S% t+T J, q ± 
being the wavevector component, perpendicular to the 
chain. In MF this is w (Jj_(0) - Jx(q±))(S} 2 /2 < 
(9/2)J±(S) 2 , where J±(q ± ) is the Fourier transform 
of the interchain coupling, and can be neglected for 
fJi(S') 2 <C — |i?Gs(l ~ cosg). Finally, the interchain 
dispersion appears, which changes the magnon average 
energy in the denominator. In the random phase ap- 
proximation (RPA), e{q,q ± ) w e(q)^J\ + x q J±(<l±) ~ 
e(q) + S(q)J±(q ± ), where e(q) and X g = 2S(q)/e(q) 
are the magnon energy and the spin susceptibility for 
a single chain. From the boundaries of the interchain 
dispersion at q — tt, measured for CSN1CI3 in Q at 
T « 10 K, and using x^.t^oaj ~ 0.74x^t=o f§], we find 
max{J±(q_ L )} — min{J±(q ± )} = 9J± s» 0.3J. Plugging 
this value in the above estimates, we conclude, that the 
purely ID expression for the static structure factor (||), 



and, by definition, the energy integrated scattering inten- 
sity, holds for CsNiCl 3 within ±5% at 0.3tt < q < 0.7tt, 
and within ±10% at 0.2tt < q < 0.8tt. 

Thus, our measurements present a detailed character- 
ization of the non-hydrodynamic part of the excitation 
spectrum in the ID S = 1 HAFM. We find that single- 
mode dispersion gradually crosses over to a narrow con- 
tinuum at q < 0.5-7T. The continuum starts below the low- 
est possible energy of the two non-interacting magnons, 
indicating their attraction. 
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